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Abstract. We present deep, high velocity resolution (~ 1.6 km sec -1 ) Giant Meterwave Radio Telescope (GMRT) 
HI 21cm synthesis images of the faint (Mb ~ —10.6) local group dwarf galaxy DDO210. We find that the HI 
distribution in the galaxy is not axi-symmetric, but shows density enhancements in the eastern and southern 
halves of the galaxy. The optical emission is lopsided with respect to the HI, most of the bright optical emission 
arises from the eastern half of the HI disk. The velocity field of the galaxy is however quite regular and shows a 
systematic large scale pattern, consistent with the rotational motion. The rotation curve for the galaxy shows a 
peak (inclination corrected) rotation velocity of only ~ 8 km sec" 1 , comparable to the random motions in the HI 
gas. After correcting for the dynamical support provided by random motions (the "asymmetric drift" correction), 
we find the corrected peak rotation velocity of ~ 16.0 km sec -1 . Mass modeling of the corrected rotation curve 
shows that the kinematics of DDO210 can be well fit with either a modified isothermal halo (with a central 
density po ~ 29 x 10~ 3 Mq pc -3 for a stellar mass to light ratio of 3.4) or an NFW halo. In the case of the NFW 
halo however, a good fit is obtained for a wide range of parameters; the halo parameters could not be uniquely 
determined from the fit. Density profiles with inner slope steeper than ~ 1.2 however provide a poor fit to the 
data. Finally, the rotation curve derived using MOND also provides a reasonable fit to the observed rotation 
curve. 
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1. Introduction 

Numerical simulations of hierarchical galaxy formation 
models such as the CDM model predict a "universal" 
cusped density core for the dark matter halos of galaxies 
(e.g. Navarro et al. 1996). A cusped density core corre- 
sponds to a steeply rising rotation curve. The observed 
kinematics of galaxies can hence be used to test such 
numerical models of galaxy formation. Dwarf low sur- 
face brightness (LSB), galaxies are best suited for such 
a test, since they, unlike larger galaxies, are known to be 
dark matter dominant. In large spiral galaxies, both gas 
and stars make significant contributions to the total mass, 
particularly in the inner regions of the galaxy. Since the 
exact contribution of the baryonic material to the total 
mass of the galaxy depends on the unknown mass to light 
ratio of the stellar disk, it is generally difficult to unam- 
biguously determine the density profile of the dark mat- 
ter halo in the central regions of large spirals. In dwarf 
LSB galaxies on the other hand, since the stellar disk is 
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generally dynamically unimportant, the central halo den- 
sity can be much better constrained. Interestingly, the ob- 
served rotation curves of dwarf galaxies generally indicate 
that their dark matter halos have constant density cores 
(e.g. Weldrake et al. 2003) unlike the cusped density cores 
predicted in numerical simulations. Another prediction of 
the numerical simulations is that the density of the dark 
matter halo is related to the background density at the 
time of the halo formation. Since the smallest galaxies 
form first in such models, these galaxies are expected to 
have the largest halo densities. The determination of the 
shapes and characteristic densities of the dark matter ha- 
los of the faintest dwarf galaxies is hence a particularly 
interesting problem. However, a major stumbling block in 
such programs is that it is currently controversial whether 
very faint dwarf irregular galaxies show systematic rota- 
tion or not. Lo et al. (1993), in a study of the kinemat- 
ics a sample of nine faint dwarfs (with Mb ~ —9.0 to 
Mb ~ —14.0) found that most of them were character- 
ized by chaotic velocity fields. However, as pointed out by 
Skill man (1996)[ Lo et al.'s observations had limited sensi- 
tivity to faint extended emission which is likely to have led 
to an underestimation of the rotation velocities. Further, 
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from a recent high sensitivity and high velocity resolu- 
tion study of the dwarf irregular galaxy Camelopardalis B 



(Cam B), Begum et al. (2003) found that inspite of being 
very faint (Mb ~ —10.8), the galaxy shows systematic ro- 
tation. Do all faint dwarf irregular galaxies have a rotating 
HI disk or is Cam B a special case? What is the dark mat- 
ter distribution in these very faint galaxies? In this paper 
we discuss these questions in the specific context of the 
local group dwarf galaxy DDO210. 

DDO210, the faintest known (M B ~ -10.6) gas 
rich dwarf galaxy in our local group, was discovered by 
van den Berg (1959) and later detected in an HI 21 cm 
survey by Fisher and Tully (1975). Fisher and Tully (1979) 
assigned a distance of 0.7 Mpc to it, based on its proxim- 
ity to NGC 6822 both on the sky and in velocity. On the 
other hand, Greggio et al. (1993) based on the colour- 



magnitude (C-M) diagram of DDO210, estimated its dis- 
tance to be 2.5 Mpc. However, recent distance estimates 
for DDO210 give distances closer to the original estimate 
of Fisher and Tully (1979). | Lee et al. (1999)1 based on the 
I magnitude of the tip of the red giant branch, estimated 
the distance to DD210 to be 950±50 kpc. This estimate is 
in excellent agreement with the value of 940±40 kpc de- 
rived recently by Karachentsev et al. (2002) using HST 



observations. At this distance, DDO210 would be a mem- 
ber of the local group. 

DDO210 is classified as a dlr/dSph or "transition 
galaxy" , with properties intermediate between dwarf ir- 
regulars and dwarf spheroidals (Mateo 1998). For exam- 
ple, in spite of containing a significant amount of neutral 
gas, DDO210 shows no signs of ongoing star formation. 
Ha imaging detected a single source of line emission in 
the galaxy; however follow up observations of this emis- 
sion suggests that it does not arise in a normal HII re- 
gion, but probably comes from dense outflowing material 
from an evolved star (van Zee et al. 1997). Consistent 
with the lack of ongoing starformation, the CM diagram 
of DDO210 shows that the brightest stars in the galaxy are 
the faintest among the brightest stars in all known dwarf 
irregular galaxies of our local group (Lee et al. 1999). 

Prior to this work, there have been two HI intcrfero- 
metric studies of DDO210, both of which used the VLA. 



This galaxy was a part of the sample of Lo et al. (1993) 



discussed above, and has also been recently re-observed 
with a high velocity resolution in the Cs array (Young et 
al. 2003). Although fYoung et al. (2003)| noted that DD210 
showed a systematic large scale velocity gradient, no at- 
tempt was made to derive a rotation curve for the galaxy 
since their study focused mainly on the local connection 
between the ISM and star formation . 

We present here deep, high velocity resolution (~ 
1.6 km/s) Giant Meterwave Radio Telescope (GMRT) ob- 
servations of the HI emission from DDO210 and use them 
to study the kinematics of the neutral gas in this galaxy. 
The rest of the paper is divided as follows. The GMRT 
observations are detailed in Sect. [5J while the results are 
presented in discussed in Sect. Throughout this paper 



we take the distance to DDO210 to be 1.0 Mpc, and hence 
its absolute magnitude to be Mb ~ —10.6. 

2. Observations 

The GMRT observations of DDO210 were conducted from 
13—15 July 2002. The setup for the observations is given 
in Tabled Absolute flux calibration was done using scans 
on the standard calibrators 3C48 and 3C286, one of which 
was observed at the start and end of each observing 
run. Phase calibration was done using the VLA calibrator 
source 2008-068 which was observed once every 30 min- 
utes. Bandpass calibration was done in the standard way 
using 3C286. 

The data were reduced using standard tasks in classic 
AIPS. For each run, bad visibility points were edited out, 
after which the data were calibrated. Calibrated data for 
all runs was combined using DBCON. The GMRT does 
not do online doppler tracking - any required doppler 
shifts have to be applied during the offline analysis. 
However since the differential doppler shift over our ob- 
serving interval is much less than the channel width, there 
was no need to apply an offline correction. 

The GMRT has a hybrid configuration (Swarup et al. 
1991) with 14 of its 30 antennas located in a central com- 
pact array with size « 1 km (ss 5 kA at 21cm) and the 
remaining antennas distributed in a roughly "Y" shaped 
configuration, giving a maximum baseline length of « 25 
km (w 120 kA at 21 cm). The baselines obtained from an- 
tennas in the central square are similar in length to those 
of the "D" array of the VLA while the baselines between 
the arm antennas are comparable in length to the "B" 
array of the VLA. A single observation with the GMRT 
hence yields information on both large and small angular 
scales. Data cubes were therefore made at various (u,v) 
ranges, including 0-5 kA, 0-10 kA, 0-15 kA and 0-20 
kA using uniform weighting. At each (u,v) range, a cir- 
cularly symmetric gaussian taper with a FWHM equal 
to 80% of the (u,v) range was applied inorder to reduce 
the sidelobes of the synthesized beam. The angular reso- 
lutions obtained for the various (u,v) ranges listed above 
were 44" x 37", 29" x 23", 20" x 15" and 12" x 11" re- 
spectively. All the data cubes were deconvolved using the 
the AIPS task IMAGR. 

The HI emission from DDO210 spanned 22 channels 
of the spectral cube. A continuum image was made using 
the average of remaining line free channels. No contin- 
uum was detected from the disk of DDO210 to a 3cr flux 
limit of 1.2 mJy/Bm (for a beam size of 30 x 23 ). We 
also checked for the presence of any compact continuum 
sources in the disk of DDO210 by making a high resolution 
(3.8 x 3.0 ) map - no sources associated with the disk of 
DDO210 were detected down to a 3a limit of 0.5 mJy/Bm. 

The line profiles were examined at various locations 
in the galaxy and were found (to zeroth order) to be 
symmetric and single peaked. Although double gaussians 
or Gauss-Hermite polynomials provide a better fit to the 
data (particularly in the higher column density regions), 
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Table 1. Parameters of the GMRT observations 



Parameters 



Value 



RA(2000) 

Declination(2000) 

Central velocity (heliocentric) 

Date of observations 

Time on source 

Total bandwidth 

Number of channels 

Channel separation 

FWHM of synthesized beam 

RMS noise per channel 



20 h 46 m 53.0 a 

-12°50'57" 

-139.0 km sec" 1 

13-15 July 2002 

16 hrs 

1.0 MHz 

128 

1.65 km sec -1 

44" x 37", 29" x 23", 

20" x 15", 12" x 11" 

2.2 mjy, 1.8 mjy, 1.6 mjy 

1.4 mjy 



the peak velocity of the profile in these regions matches 
(within the errorbars) with the intensity weighted mean 
velocity. Since we are interested here mainly in the sys- 
tematic velocities, moment maps provide an adequate de- 
scription of the data. Moment maps were therefore made 
from the data cubes using the MPS task MOMNT. To 
obtain the moment maps, lines of sight with a low signal 
to noise ratio were excluded by applying a cutoff at the 
3cr level, (cr being the rms noise level in a line free chan- 
nel), after smoothing in velocity (using boxcar smoothing 
three channels wide) and position (using a gaussian with 
FWHM ~ 2 times that of the synthesized beam) . Maps of 
the velocity field and velocity dispersion were also made in 
GIPSY using single Gaussian fits to the individual profiles. 
The velocity field produced by gaussian fitting is in rea- 
sonable agreement with that obtained from moment anal- 
ysis. However, the moment2 map systematically underes- 
timates the velocity dispersion. This can be understood 
as the effect of the thresholding algorithm used by the 
MOMNT task to identify the regions with signal. The ve- 
locity dispersion from gaussian fitting to the lines profiles 
was found to be ~ 6.5 km sec -1 and showed no measurable 
variation across the galaxy. This lack of substantial varia- 
tion of o across DDO210 is typical of dwarf irregular galax- 
ies, although the observed velocity dispersion for DDO210 
is lower than the typical value observed in dwarf irregular 
Lake et al. (1990)| [Skillman et al. (1988) I 



galaxies (e.g. 

A lower value of HI dispersion for DDO210 was also noted 



by Young et al. (2003) 



3. Results and Discussion 

3.1. HI distribution 

The global HI emission profile of DDO210, obtained from 
44" x 37" data cube, is shown in Fig. ^ A Gaussian 
fit to the profile gives a central velocity (heliocentric) 
of -139.5 ± 2.0 km s" 1 . The integrated flux is 12.1 ± 
1.2 Jy km s . These are in good agreement with the val- 
ues of -140±2.0 km s _1 and 11.5± 1.2 Jy km s _1 obtained 
from single dish observations (Huchtmeier & Richter 
1986). A good agreement between the GMRT flux and 



the single dish flux shows that no flux was missed because 
of the missing short spacings in the interferometric obser- 
vations. The velocity width at 50 % level of peak emission 
(AV50) is found to be 19.1±1 km s _1 , which again is in ex- 
cellent agreement with the AV50 value of 19.0 km s _1 from 
the single dish observations. The HI mass obtained from 
the integrated profile (taking the distance to the galaxy to 
be 1.0 Mpc) is 2.8 ± 0.3 x 10 6 M Q , and the M m /L B ratio 
is found to be ~ 1.0 in solar units. 

Fig.[3shows the integrated HI emission from DDO210 
at 44" x 37" resolution, overlayed on the digitised sky 
survey (DSS) image. The HI isodensity contours are elon- 
gated in eastern and southern half of the galaxy i.e. the 
density is enhanced in these directions. These density en- 
hancements are highlighted in Fig. 01 which shows the in- 
tegrated HI emission at high resolution (12" x 11"). Note 
that the faint extended emission seen in the low resolution 
map is resolved out in this image. 

The optical emission in the galaxy shows two com- 
ponents, a central bright compact component and an 
outer faint extended component; both elongated in the 
east- west direction (Lee et al. 1999). However, the cen- 
ters for the two components do not overlap; the bright 
component lies more eastwards than the center of ex- 
tended component. From a deep C-M diagram of DDO210, 
|Tolstoy et al. (2000)| found a n umber of faint s tars older 
than several Gyr in the galaxy, van Zee (2000) estimated 
an inclination of 68 degrees for DDO210 from her obser- 
vations. Lee et al. (1999)[ based on deeper observations, 
estimated the inclination of the faint extended component 
to be ~60 degrees, but noted that the presence of several 
bright foreground stars makes this estimate of inclination 
very uncertain. 

From a visual inspection of the overlay (Fig. |5J), the 
eastern HI density enhancements seem to correlate with 
the optical emission in the galaxy. As can be seen the 
bright component of the optical emission is not centered 
on the HI emission, but is instead offset to the east. On 
the other hand, no optical emission is seen along the 
southern HI density enhancement. As a further check for 
optical emission associated with the southern HI den- 
sity enhancement, we computed number counts of stars 
in that region using t he deep optical image obtained 
by Tolstoy et al. (2000) No significant increase in the star 
counts from the background was detected. 

As seen in the Fig. the faint extended emission in 
the lower resolution HI distribution is less distorted by the 
presence of the HI density enhancements. Hence, an esti- 
mate of the morphological center, position angle (PA) and 
inclination of the galaxy was obtained by fitting elliptical 
annuli to the 44" x 37" resolution HI distribution. The 
estimated HI morphological center lies close to the cen- 
ter of the faint extended optical component (but as noted 
above is offset from the center of the bright optical emis- 
sion). The variation of the fitted morphological PA with 
the galactocentric radius is given in Fig. [5] The inclina- 
tion estimated from the outer HI contours (assuming the 
intrinsic shape of the HI disk to be circular) is 27±7.0 de- 
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grees. The inclination in the inner regions of the galaxy 
could not be constrained reliably due to the significant 
distortion in the HI contours. For the same reason, ellipse 
fitting to the higher resolution HI data (where the smooth 
emission is resolved out) is not reliable. 
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Fig. 1. HI profile for DDO210 obtained from 44 x 
37 data cube. The channel separation is 1.65 km s _1 . 
Integration of profile gives a flux integral of 12.1 Jy 
km s _1 and an HI mass of 2.8 x 1O 6 M . The dashed line 
shows a gaussian fit to the profile. 

The inclination derived from the HI distribution is 
significantly different from the optical inclination of the 
galaxy. It is likely that the optical emission does not arise 
from an axi-symmetric disk, but instead that the stars are 
concentrated in a patchy elongated region in DDO210. 
Hence, in the absence of any other reliable estimate, the 
value of inclination obtained from the outer HI contours 
was assigned to the whole galaxy. The deprojected HI ra- 
dial surface density profiles were then obtained by aver- 
aging over elliptical annuli in the plane of the galaxy. The 
profiles derived from the 44" x 37" and 29" x 23" resolu- 
tion HI distributions are given in Fig. U3 As can be seen, 
the two distributions are in reasonable agreement. The 
flux integral estimated from the 29" x 23" resolution HI 
moment map was found to be ~11% smaller than that 
estimated from the 44" x 37" resolution HI distribution. 
For the rest of the analysis, we will use the 44" x 37" HI 
profile. This profile is well represented by a gaussian, i.e. 
we have 



£ ffl (r) = So x ' 2r ° 

with S =13.9±1.0 M Q pc- 2 , 
r = 46" (~ 0.22kpc). 



c=-4.0'' 



o.o2k pc ; 



(i) 

and 
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Fig. 2. The optical DSS image of DDO210 (greyscales) 
with the GMRT 44" x 37" resolution integrated HI emis- 
sion map (contours) overlayed. The contour levels are 0.01, 
0.151, 0.293, 0.434, 0.576, 0.717, 0.854, 0.999, 1.141, 1.283, 
1.424, 1.565, 1.797, 1.848 Jy/Bm km s^ 1 
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Fig. 3. The HI surface density profile derived from the HI 
distribution at 29 x 23 (squares) and 44 x 37 (trian- 
gles) resolution. A gaussian fit to 44 x 37 HI distribution 
is shown superimposed. 
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Fig. 4. Integrated HI emission map of DDO210 (grey 
scale and contours) at 12 x 11 resolution. The contour 
levels are 0.001, 0.154, 0.293, 0.434, 0.576, 0.717, 0.858, 
0.999, 1.141 Jy/Bm km s" 1 
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Fig. 5. The variation of position angle (PA) with the 
galaetocentric distance. Circles represent the morphologi- 
cal PA obtained from ellipse fit to the HI distribution at 
44 x 37 resolution. Triangles represent the kinematical 
PA derived from 44 x 37 resolution velocity field. 

3.2. HI Kinematics 

The velocity field of DDO210 derived from the moment 
analysis of 29" x 23" resolution data cube is shown in 
Fig. The velocity field is regular and a systematic ve- 
locity gradient is seen across the galaxy. Our velocity 
field differs significantly from the velocity field derived by 



Fig. 6. The HI velocity field of DDO210 at 29 x 23 
resolution. The contours are in steps of 1 km sec -1 and 
range from —145.0 km sec" 1 to —133.0 km sec -1 . 



ity field is, to zeroth order, consistent with that expected 
from rotation. On the other hand, the velocity field derived 
by Lo et al. (1993) (based on a coarser velocity resolution 
of ~ 6 km s _i ) is chaotic. This difference in the observed 
kinematics suggests that high velocity resolution and high 
sensitivity is crucial in determining the systematic gradi- 
ents in the velocity field of faint galaxies like DDO210. 
A similar systematic velocity pattern was also found in 
the recent high velocity resolution study of DDO210 by 
Young et al. (2003) 



Lo et al. (1993) The systematic pattern seen in our veloc- 



however, those authors concentrated 
on trying to determine the physical conditions of the emit- 
ting gas and not the large scale kinematics of the galaxy. 

The velocity field in Fig. H3 also shows signa- 
tures of both warping and kinematical lopsidedness 
IjSwaters et al. 1999|l . i.e. the kinematical major axis is 
not a straight line but shows twists and the isovelocity 
contours in the eastern half of the galaxy are more closed 
than those in the western half. Both of these kinemat- 
ical peculiarities become more prominent in the higher 
spatial resolution velocity field (see Fig. [SJ. While the 
density enhancements in the HI also become more promi- 
nent at higher spatial resolutions, there does not seem to 
be any particular correlation between the HI density en- 
hancements and these kinematical peculiarities. 

3.3. HI Rotation Curve 

Given the lack of correlation between the density enhance- 
ments in the HI distribution and the global kinematical 
peculiarities in DDO210, it would be reasonable to assume 
that these density enhancements follow the same kinemat- 
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Fig. 7. The rotation curves derived from the intensity 
weighted velocity held at various resolutions. Crosses, tri- 
angles, circles and squares show the rotation velocity de- 
rived from the 12"xll",20"xl5", 29"x23" and 44"x37" 
resolution respectively. The adopted hybrid rotation curve 
is shown by a solid line. 

ics as that of the more diffuse extended emission. Under 
the further assumption that the kinematical peculiarities 
noted above are not important (to the zeroth order; see 
also Sect. 13. 5H . rotation curves were derived by fitting the 
usual tilted ring model to the HI velocity fields at various 
resolutions (including 44" x 37", 29" x 23", 20" x 15" and 
12" x 11") using the GIPSY task ROTCUR. 

First, the kinematical center and the systemic velocity 
(V sys ) of DDO210 were obtained from a global fit to the 
velocity fields at various resolutions. The values of V sys 
derived from the velocity fields matched within the er- 
rorbars and also matched with the value obtained from 
a gaussian fit to the global HI profile. The kinematical 
center derived from a global fit to the velocity fields at 
44" x 37" and 29" x 23" resolution matched within the 
errorbars and also matched with the morphological cen- 
ter derived from the 44" x 37" resolution HI distribution. 
However, because of distorted morphology of DDO210 at 
high spatial resolutions, an attempt to derive the kine- 
matical center from a global fit to the velocity fields at 
these resolutions did not yield reliable results. Hence, the 
center for the higher resolution velocity fields was fixed to 
a value obtained from the lower resolution velocity fields. 

Keeping V sys and the kinematical center fixed to the 
values obtained from the global fit, the position angle (PA) 
of the galaxy was derived, using tilted ring model, by 
breaking up the galaxy into elliptical annuli (each of width 
half that of the synthesized beam). The variation of the 
derived kinematical PA with the galactocentric radius for 
the 44" x 37" resolution velocity field is given in the Fig. [5] 
The value of kinematical PA at all galactocentric radii 
matched with the morphological PA estimated from the 
44" x 37" resolution HI distribution. Attempts to derive 



the kinematical position angle at the higher spatial resolu- 
tions did not yield reliable estimates. Similarly, attempts 
to derive the kinematical inclination of DDO210 did not 
give reliable results at any spatial resolution. Hence, in the 
absence of any other reliable estimate, the kinematical in- 
clination of the HI disk was fixed to the value obtained 
from the HI morphology. Finally, the rotation curve of 
the galaxy was computed keeping all parameters, except 
the circular velocity V c , in each elliptical annuli fixed. 

Fig. [3 shows the rotation curves derived from the dif- 
ferent spatial resolution velocity fields. As can be seen, the 
20" x 15" and 12" x 11" resolution velocity fields lead to the 
rotation curves that are significantly steeper in the inner 
regions of the galaxy than those computed from the lower 
resolution 44" x 37" and 29" x 23" velocity fields; this could 
be a result of beam smearing. On the other hand, the rota- 
tion curves from each resolution agree in the outer regions 
of the galaxy. This suggests that the effect of beam smear- 
ing is significant only in the inner regions of the galaxy, 
where the rotation velocity is increasing with the galacto- 
centric radius. Further, the rotation curves derived from 
the 12" x 11" resolution velocity field agrees with that ob- 
tained from the 20" x 15" resolution data suggesting that 
beam smearing effects are no longer significant at such 
high resolutions (which corresponding to a linear scale of 
~ 70pc). The final rotation curve that we adopt for the 
rest of the analysis is shown as a solid line in Fig. [3 The 
rotation velocities derived from the 20" x 15" resolution 
velocity field are used in the inner regions of the galaxy 
(upto ~ 50"), while, the outer points are taken from the 
29" x 23" and the 44" x 37" resolution data. 

As discussed above, the rotation curves have been de- 
rived under the assumption that the HI density enhance- 
ments follow the same kinematics as that of the extended 
emission. The validity of this assumption could be tested 
by checking whether the derived rotation curves can repro- 
duce the observed kinematics in the galaxy. Hence, model 
velocity fields were made at 44" x 37", 29" x 23" and 
20" x 15" resolutions using the rotation curves and other 
kinematical parameters, derived from each resolution with 
a tilted ring model, using the task VELFI in GIPSY. The 
model velocity field at 20" x 15" resolution is shown in 
Fig- [HI As can be seen from the figure, the model provides 
a reasonable match to the observed velocity field. Note 
however that the model fails to reproduce the kinematical 
lopsidedness and the twist in the kinematical major axis. 
This is to be expected, since neither of these two features 
was incorporated in the model. A similar match between 
the observed and the model velocity fields was also found 
for the two lower resolution velocity fields. 

The maximum rotation velocity (inclination corrected) 
of DDO210 is ~ 8.0 km s" 1 , which is comparable to the 
velocity dispersion observed in the HI gas, i.e. random mo- 
tion in the gas provides significant dynamical support to 
the HI disk. Since the pressure support contributes signifi- 
cantly to the dynamics of DDO210, the observed rotation 
velocities underestimate the total dynamical mass. The 
rotation velocity hence has to be corrected for the pres- 
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Fig. 8. [A] The observed velocity field of DDO210 at 20" x 15" arcsec resolution. [B]The model velocity field derived 
from the rotation curve at 20" x 15" resolution. The contours are in steps of 1 km sec -1 and range from -144.0 km 
sec -1 to -134.0 km sec -1 . 



sure support before one can derive a mass model for the 
galaxy; this correction (generally called the "asymmetric 
drift" correction) is given by: 

v 2 c =v 2 ~rxa 2 (In Shi) + ^(m^x 2 ) - £(ln2fc,)] , (2) 



dr 



dr 



where v c is the corrected circular velocity, v a is the ob- 
served rotation velocity, a is the velocity dispersion, and 
h z is the scale height of the disk. Strictly speaking, asym- 
metric drift corrections are applicable to collisionless stel- 
lar systems for which the magnitude of the random mo- 
tions is much smaller than that of the rotation velocity. 
However, it is often used even for gaseous disks, where the 
assumption being made is that the pressure support can 
be approximated as the gas density times the square of 
the random velocity. In the absence of any measurement 
for h z for DDO210, we have assumed d(ln(h z ))/dr = 
(i.e. that the scale height does not change with radius). 
Also, using the fact that a is constant across the galaxy, 
we get: 



r xa 2 [f (In Shi)]. 
L dr J 



(3) 



Using the fitted Gaussian profile to the radial surface den- 
sity distribution, (see eqn^) we obtain 



vl + r(r - c) x cr 2 /r 2 . 



(4) 



The observed HI velocity dispersion was corrected for 
the instrumental broadening as well as for the broaden- 
ing due to the velocity gradient over the finite size of the 
beam. The applied correction is 

a t 2 rU e = ^ 2 bs-A^ 2 -^ 2 (V Uo ) 2 , 

where CTtmc is the true velocity dispersion, Av is the 
channel width, b characterizes the beam width (i.e. the 



beam is assumed to be of form e~ x l h ) and v a is the ob- 
served rotation velocity. After putting the appropriate val- 
ues in the above equation, we get of rue ~ 36.0 km 2 sec -2 . 
Finally, substituting this value back into the Eqn. J3J the 
"asymmetric drift" corrected curve obtained is given as a 
dotted line in Fig. |5J 

As discussed above, the "asymmetric drift" correction 
was derived under the assumption that the disk scale 
height of DDO210 does not change with radius. To quan- 
tify the effect of this assumption, the correction was re- 
calculated assuming a linear increase of 100% in the scale 
height, from q =0.25 at the center to q =0.5 at the edge 
of the galaxy. The change in the asymmetric drift cor- 
rection was found to be < 1%. The assumption that the 
scale height is constant hence have a negligible effect on 
the derived halo parameters. 

3.4. Mass Model 

In this section we use the "asymmetric drift" corrected 
rotation curve derived in the last section to derive mass 
models for DDO210. 

As discussed in Sect. 13.11 the optical emission is not 
axi-symmetric but is instead patchy and elongated. For 
the purpose of mass modeling however, we approximate 
the optical emission to be axi-symmetric, but we return to 
this issue in Sect 13.51 Assuming that the optical emission 
was axi-symmetric, Lee et al. (1999) estimate the B band 
scale length (a B ) of DDO210 to be 35.0" (~ 0.17 kpc). 
This estimate of olb agrees to within 10% with the value 
of 39" estimated by|van Zee (2000)|| Lee et al. (1999)|also 



found the colours of the galaxy to be nearly constant 
with the galactocentric radius. Hence, we computed the 
contribution of the stellar mass to the observed rotation 
curve by assuming it to be an exponential disk of con- 
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Fig. 9. The hybrid rotation curve (dashes) and the rota- 
tion curve after applying the asymmetric drift correction 
(dots). 



stant M/Lb ratio (Tb) and with an intrinsic thickness 
ratio (qo) of 0.25. For the vertical density distribution of 
the stellar disk, we assumed a sech 2 (z/,zo) profile, with Zq 
independent of galactocentric radius. These are reasonable 
assumptions for disk galaxies (e.g. van der Kruit & Searle 
1981, de Grijs & Peletier 1997). In the absence of any prior 
knowledge of , it was taken as a free parameter in the 
mass modelling. 

The contribution of the HI mass to the observed rota- 
tion velocities was calculated using the HI surface density 
profile derived from 44" x 37" HI distribution. To correct 
for the mass fraction of Helium, the HI mass was scaled 
by a factor of 1.25. A search for molecular gas in DDO210 
gave negative result (Taylor et al. 1998), hence, the contri- 
bution of molecular gas to the rotation curve was ignored. 
We also neglected the contribution of ionized gas, if any. 
Not much is known about the vertical distribution of gas in 
dwarf irregular galaxies, however there is some evidence of 
similar vertical density distribution of HI and stellar disk 
(Bottema et al. 1986). Hence, for the HI disk we again 
assumed qo of 0.25, with a vertical density distribution 
profile of sech 2 (z/zo)- The rotation velocities for the HI 
and the stellar components were then computed using the 



formulae given by Casertano (1983) 



For the dark matter halo we considered two types of 
density profiles, viz. the modified isothermal profile and 
the NFW profile. The modified isothermal density profile 
is given by: 

PisoM = p [l + (r/rcfr 1 , 

where, po is the central density of the halo and r c is the 
core radius. The corresponding circular velocity is given 
by: 



Kr) = v /4 7 rGp r 2 [l-^ta n - 1 (^)]. 
The NFW halo density profile is given by: 
PnfwW = Pi/[(r/r s )(l + r/r s ) 2 ], 



where, r s is the characteristic radius of the halo and 
pi is the characteristic density. The circular velocity can 
be written as: 

/ \ / lnfl+ca;) — cx/(l+cx) 

V{X) = V 200 ^ a V Nl+ ;)-c/(l+c)] ' 

where, c = rzoo/r s , x = r/r^oo'i T2oo is the distance at 
which the mean density of the halo is equal to 200 times 
the critical density and t>2oo is the circular velocity at this 
radius. 
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Fig. 10. Mass models for DDO210 using the corrected 
rotation curve. The points are the observed data. The to- 
tal mass of gaseous disk (dashed line) is 3.6 x 10 6 MQ.Thc 
stellar disk (short dash dot line) has Tb = 3.4, giving 
a stellar mass of 9.2 x 10 6 Af©. The best fit total rota- 
tion curve for the constant density halo model is shown 
as a solid line, while the contribution of the halo itself is 
shown as a long dash dot line (the halo density is den- 
sity po = 29 x 10 _3 M Q pc -3 ). The best fit total rotation 
curve for an NFW type halo, using Tb = 0.5, c=5.0 and 
^200 =38.0 km s~ris shown as a dotted line. See the text 
for more details. 

First, we consider a modified isothermal dark halo 
mass model. Since the rotation curve keeps on rising till 
the last measured point (see Fig. EJl , the core radius r c 
could not be constrained. This essentially reduces a mod- 
ified isothermal dark halo fit into a constant density halo 
fit. The remaining free parameters are the central density 
of the halo, po, and the mass to light ratio of the stellar 
disk, Tb- 

Fig. shows the best fit mass model for a con- 
stant density halo. The best fit gave Tb of 3.4±0.5 and 
po = 29 ± 5 x 10" 3 Mq pc~ 3 . The observed B-V for 
DDO210 is -0.25 (Lee et al. 1999), which, (from the low 
metallicity Bruzual & Chariot SPS model using a modi- 
fied Salpeter IMF, Bell & de Jong 2001), corresponds to a 
Tb of — 0.5 . If we adopt this value of Tb then the best 
fit model has po = 59 ±7 x 10~ 3 Af© pc~ 3 . However, with 
this Tb, the fit to the observed rotation curve is much 
poorer and the model curve substantially underestimates 
the observed rotation velocities at small radii. The best 
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fit model (Yfl = 3.4) gives the mass of stellar disk to be 
M*=9.2 x 10 6 Af Q .The mass of the gas disk in DDO210 is 
M gas =3.5 x 10 6 Af Q . From the last measured point of the 
observed rotation curve, we get a total dynamical mass 
of Mx=3.4x 10 7 Af Q , i.e. at the last measured point more 
than 63% of the mass of DDO210 is dark. 

A similar procedure was also tried using a dark matter 
halo of the NFW type. Keeping as a free parameter 
in the fit gave unphysical results, hence, it was kept fixed 
to a more likely value (viz. = 0.5). We found that an 
NFW halo provides a good fit to the data, for a wide range 
of values of V200 an d c. The range of parameters which 
provide acceptable fits are (v 2 oo ~20 km s -1 , c ~ 10) 
to (i>2oo ^500 km s _1 , c ~ 0.001). As an illustration we 
show in Fig. ^| the best fit rotation curve for an NFW 
halo, using = 0.5, c = 5, U200 = 38 km s -1 . As can 
be seen clearly, the NFW halo also provides a good fit to 
the data. However, the best fit values of the concentration 
parameter c, at any given V200 was found to be consistently 
smaller than the value predicted by numerical simulations 
for the ACDM universe (Marchesini et al. 2002). 

As seen above, both isothermal and NFW halo pro- 
vide a good fit to the observed kinematics of DDO210. 
Profiles steeper than NFW have also been proposed by 
some N-body simulations (e.g. Moore et al. 1999). In or- 
der to check whether such steep profiles are also consistent 
with the data, mass models were fit using a broader family 
of density profiles, viz. 



p(r) 



(r/r ) a [l + (r/r o yf~ a)h 



(5) 



The circular velocity corresponding to the above den- 
sity profile (see Kravtsov et al. 1998) is: 



V(r) - V t 



(r/rt) E 



[l + (r/r t ) a ] (g+b)/a 



(G) 



where r t and Vt are the effective "turnover" radius 
and velocity. The parameter "g" is related to the inner 
slope of the density profile, a by g= 1 + a/2.0, "b" is the 
outer logarithmic slope of the rotation curve while "a" 
determines the sharpness of turnover. Fits to the rotation 
curve with all three parameters left free did not converge. 
Hence, since we are primarily concerned with the slope in 
the inner regions, we fixed the parameters b and a to the 
values of 0.34 and 1.5 respectively, (which are the typi- 
cal values found for the rotation curves of dwarf galaxies; 
Kravtsov et al. 1998), while r t and Vt were left as free pa- 
rameters. Tb was fixed to a value of 0.5, as suggested by 
the observed colours in the galaxy; this also allows a mean- 
ingful comparison of derived mass models for a family of 
density profiles. We found that the reduced x 2 for the fit 
continuously increases as the profile gets steeper. For com- 
parison, fixing g to 0.5 (corresponding to a =1.0; NFW 
profile) gave reduced \ 2 =0.4 while g of 0.4 (corresponding 
to a =1.2) gave reduced % 2 =0.7. At the extreme, fixing 
g to 0.2 (a =1.6) (which substantially over-predicts the 



observed rotation velocity at small radii, while underesti- 
mates the velocity at large radii) gave reduced x 2 =2.5. 
Note that since there are heuristics involved in computing 
the error bars on v c , it is not possible to rigorously trans- 
late the minimum x 2 value into a confidence interval for 
the parameters of the fit. However, a lower x 2 value does 
imply a better fit (see also the discussion in van den Bosh 
& Swaters 2001). Fig. ITT1 shows the best fit mass model 
with a = 1.2. 
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Fig. 11. The best fit MOND rotation curve (solid line) 
to the "asymmetric drift" corrected rotation velocities 
(points). Also shown in the figure is the best fit total rota- 
tion curve for a dark halo with a = 1.2 inner slope (dotted 
line). See the text for more details. 



So far we have assumed that the discrepancy between 
the dynamical mass, estimated from the "asymmetric 
drift" corrected rotation curve, and the luminous mass 
in DDO210 can be explained by considering an extended 
dark halo around the galaxy. An alternate explanation 
for this discrepancy is that the dynamics becomes non- 
Newtonian in the limit of low acceleration, i.e. the MOND 
theory (Milgrom 1983). We have also tried fitting the ro- 
tation curve using the MOND prescription. In this fit, Tb 
and an (the critical acceleration parameter) were taken 
as free parameters. Fig. 1111 shows the best fit MOND ro- 
tation curve. As can be seen, the MOND rotation curve 
agrees well with the observed curve in the inner regions 
of the galaxy, while it underestimate the observed curve 
(by up to 2.0 kms -1 ) in the outer regions. The best fit 
model gave Tg of 0.4±0.2 and ao of 1.7±0.3 (in units 
of 10~ 8 cms -2 ) with reduced x 2 =0.7. The best fit value 
of Ts agrees with the value expected from the observed 
colours in DDO210. Also, the best fit value of ao is con- 
sistent with the mean value of ao found for other, brighter 
galaxies (Kent (1987); see also Milgrom 1988). DDO210 
is the faintest known dwarf irregular galaxy for which the 
MOND prescription provides a reasonably good fit to the 
observed kinematics. 
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Recall that all the above mass models have been de- 
rived by assigning an inclination of 30 degrees (obtained 
from the outer HI contours, as discuss in Sect l3.1|) . Inorder 
to estimate the effect of using an erroneous inclination on 
the derived halo parameters, mass modeling was also tried 
with two other values of inclination viz. 60 degrees (i.e. 
that estimated from the optical isophotes) and 45 degrees 
(a value between the optical and HI inclination). 

The best fit mass model for a constant density halo, 
using an inclination of 60 degrees, gave Ts of 1.4±0.2 and 
po = 27 ± 2.0 x 10" 3 Af Q pc" 3 with reduced x 2 = 0.3. On 
the other hand, an inclination of 45 degrees gave the best 
fit model with T B of 1.8±0.2 and p Q = 31 ±2.0 x 10~ 3 M Q 
pc -3 with reduced x 2 = 0.2. Recall that an inclination of 
30 degrees gave the best fit with of 3.4±0.5 and po = 
29 ± 5.0 x 1O~ 3 M pc -3 with reduced x 2 = 0-4. As can 
be seen, the central halo density is relatively insensitive to 
the assumed inclination. However, the best fit Tb changes 
significantly with the assumed inclination. 

For the NFW halo fit, keeping as a free parameter 
in the fit gave unphysical results, hence, it was kept fixed 
to a value of 0.5, as suggested by the observed colours in 
the galaxy. The best fit NFW model with an inclination of 
60 degrees gave a reduced x 2 = 4.0, while for 45 degrees 
a reduced % 2 = 1.6 was obtained from the best fit. As 
can be seen, the NFW halo provides a poor fit to the 
data for higher values of inclination (for comparison, an 
inclination of 30 degrees gave the best fit NFW model with 
a reduced % 2 — 0.5). Hence, an NFW halo can be ruled 
out for DDO210, if the inclination of the galaxy is higher 
than 30 degrees. 

3.5. Discussion 

As discussed in the Sect. 13.21 DDO210 shows signatures of 
twisting of the kinematical major axis as well as of kine- 
matical lopsidedness. This may be related to the lopsided- 
ness seen in the optical disk; the bulk of the optical emis- 
sion (i.e. the "bright component" discussed in Sect. 13 - lfl 
lies in the eastern half of the HI disk. A similar pattern 
of kinematical lopsidedness and lopsidedness in the opti- 
cal disk is also seen in other galaxies (e.g. M101, Bosma 
et al. 1981). To quantify the kinematical lopsidedness in 
the galaxy, rotation curves were derived at various spatial 
resolutions separately for the approaching and receding 
halves of the galaxy. Fig. El shows the hybrid rotation 
curves for the approaching and receding sides as well as 
for the galaxy as a whole. The rotation curves obtained 
from the two lowest resolution velocity fields (44 x 37 
and 29 x 23 ) matched within the errorbars with the one 
derived from the whole galaxy. On the other hand, for the 
higher resolution data, the rotation curves were found to 
be different (at ~ 2.0 km s _1 level) from the rotation curve 
derived from the whole galaxy. While part of this differ- 
ence may be due to the kinematical lopsidedness, it is also 
possible that part of it arises as a result of the difference in 
the sampling of the velocity field between the two sides of 



the galaxy. Because of the distorted morphology of the HI 
distribution, the distribution of data points available to 
estimate rotation velocity are different for the approach- 
ing and receding halves. This effect gets more important 
at the higher spatial resolutions and could also contribute 
to the observed differences in the rotation curves. Further, 
as can be seen in Fig. ^| the difference between the rota- 
tion curve derived using the whole galaxy and the rota- 
tion curves derived separately for the two halves are small 
compared to both the error bars on the average curve as 
well as the magnitude of the asymmetric drift correction. 
Ignoring the kinematical lopsidedness for the mass mod- 
eling is hence probably reasonable. 
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Fig. 12. Filled triangles show the rotation velocities de- 
rived from the whole galaxy while circles and squares rep- 
resent the rotation velocities derived from the receding 
and the approaching side respectively. 

The derived rotation curve can also be used to put 
limits on the dynamical ages of the density enhancements 
in the HI distribution. The derived rotation curve for the 
galaxy (see Fig.[7J) is flat from a galactocentric radii of 0.2 
kpc to 0.5 kpc. The time scales required for one rotation 
in the inner regions of the galaxy is ~ 160 Myr, while the 
rotation period at the edge of the galaxy is ~ 400 Myr. 
Hence, this differential rotation will wind up these density 
enhancements on a time scale of a few 100 Myrs. These 
density enhancements are hence likely to be recently de- 
veloped perturbations in the HI disk. The origin of these 
perturbations is unknown, and it also surprising that de- 
spite their presence and relative dynamical youth there is 
not much evidence for star formation in the galaxy. 

To conclude, we have presented deep, high velocity res- 
olution (~ 1.6 km sec -1 ) GMRT HI 21cm synthesis im- 
ages for the faint (Mb 10.6) dwarf galaxy DDO210. 

We find that the HI distribution in the galaxy is not axi- 
symmetric, but shows density enhancements in the eastern 
and southern parts of the galaxy. The velocity field of the 
galaxy is however regular, and shows a systematic large 
scale pattern, consistent with rotational motion. The high 
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velocity resolution (~ 1.6 km sec -1 ) and sensitivity of our 
observations were crucial for determining the true velocity 
field of the galaxy. The inclination angle for DDO210 is 
poorly known. From the observed velocity field, we derive 
a rotation curve for DDO210 by assigning an inclination 
derived from outer HI contours to the whole galaxy. The 
derived peak rotational velocity is comparable to the ran- 
dom motions in the HI gas. After correcting for the dy- 
namical support provided by the random motion, we find 
that the rotation curve of DDO210 can be well fit with 
either a modified isothermal halo (with a central density 
po ~ 29 x 1(T 3 M Q pc~ 3 ) or an NFW halo. We find that 
for a constant density halo, the central halo density is 
relatively insensitive to the assumed value of inclination. 
However, the best fit Tg changes significantly with the as- 
sumed inclination. On the other hand, an NFW halo can 
be ruled out for DDO210, if the inclination of the galaxy 
is higher than ~ 30 degrees. Mass models are also derived 
for a family of density profiles steeper than NFW; we find 
that density profiles with inner slopes steeper than ^1.2 
provide a poor fit to the data. 

As an alternative to the dark matter hypothesis, we 
have derived MOND fits to the rotation curve. We find 
that the rotation curve predicted by MOND provides an 
acceptable fit to the observed rotation curve. Further, 
the value of the critical acceleration ao of 1.7±0.3 xlO -8 
cms -2 is consistent with the earlier determinations of the 
value of this parameter. 
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